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INTRODUCTION 


There  is  considerable  interest  in  determining  the  electrical  properties 
of  turbulent  afterburning  exhaust  plumes  of  solid  propellant  rocket  motors. 
Modeling  calculations  have  been  employed  to  calculate  the  radar  interference 
effects  that  arise  when  communication  with  rockets  involves  transmission 
through  these  plumes.  In  a recent  paper,  Jensen  and  Pergament  have  shown 
that  for  systems  where  potassium  is  a major  rocket  exhaust  impurity,  above- 
equilibrium electron  concentrations  are  principally  formed  by  the  chemi- 
ionization  reaction 

K + Cl  - K+  + Cl"  (1) 

followed  by  the  electron  detachment  reaction 

Cl"  + H - HC1  + e"  . (2) 

However,  one  of  the  major  uncertainties  in  the  above  modeling  study  is  the 
ionization  reaction  rate  for  K + Cl,  which  is  known  only  to  an  accuracy  of 
one  or  two  orders  of  magnitude.  There  is  a similar  interest  in  the  Na  + Cl 
ionization  reaction. 

Some  information  on  both  forward  and  reverse  reaction  rates  for  alkali- 
halogen  atom  pairs  has  previously  been  obtained  in  bulk  experimental 
arrangements,  but  the  results  were  not  very  reliable,  and  no  theoretical 
work  had  been  done  on  these  species.  Improvements  in  both  the  experimental 
measurements  and  in  the  theoretical  picture  were  desirable.  An  unusual 
combination  of  experimental  and  theoretical  capabilities  to  treat  this 
problem  exist  in  the  Molecular  Physics  Center  at  SRI. 

The  particular  reactions  of  importance  for  rocket  plume  processes  have 
never  been  studied  by  beam  techniques,  although  certain  examples  of  the 
reverse  reaction  have  been  studied  at  the  FOM  Laboratory  in  Amsterdam. 

Our  unique  apparatus  for  studying  ion-ion  neutralization  has  previously  made 
it  possible  for  us  to  conduct  an  extensive  set  of  measurements  on  the 


neutralization  reactions  of  a considerable  number  of  ion  pairs  of  species 

+ - 

containing  H,  0 and  N.  Cross  sections  for  fifteen  ion  pairs  such  as  H + II 

+ - + - + “ 4"  “ + - -+■  -f  - 

N + 0 , 0 + 0 , NO  + 0 , C>2  + 0 , + °2  » °2  + °29  and  N0  + N°2 

have  been  measured  and  rate  constants  determined.  The  work  reported  here 

was  a small  effort  in  which  it  was  originally  intended  to  use  our  existing 

beams  equipment  to  extend  our  measurements  of  ion-ion  neutralization  cross 

sections  to  include  alkali  and  halogen  ion  pairs.  It  was  also  intended  to 

supplement  this  work  by  theoretical  calculations,  in  order  that  the 

achievement  of  reasonable  agreement  between  experiment  and  theory  in  the 

case  of  one  or  more  alkali-halogen  combinations  would  provide  further 

confirmation  of  both  approaches. 

Unexpected  experimental  difficulties  encountered  in  the  course  of  the 
work  have  prevented  a successful  measurement  within  the  limited  scope  of 
effort  and  funds  originally  allocated  to  this  project.  The  difficulties 
have  been  identified,  as  have  the  modifications  of  the  equipment  and  method 
that  would  be  required  for  successful  measurement. 

The  theoretical  calculation  has  been  successfully  carried  out  employing 
significant  improvements  in  computational  technique  and  in  parameters 
employed  as  compared  to  previous  theoretical  treatments  of  these  reactions. 
In  particular,  close  coupling  calculations  have  been  used  instead  of  the 
simplified  Landau-Zener  formula.  For  the  case  of  K + Cl,  our  current 
results  differ  from  the  Landau-Zener  formula  by  several  orders  of  magnitude. 

EXPERIMENTAL 

The  experimental  program  was  intended  to  take  advantage  of  existing 

facilities  with  minimal  modifications  compatible  with  a small  scale  effort. 

For  the  initial  neutralization  measurements,  we  chose  to  work  on  the  Na+  + F 

reaction  because  both  atoms  have  only  one  stable  isotope.  An  attempt  to 

35  37 

use  isotopically  mixed  beams  such  as  Cl“  + Cl  would  have  led  to 
experimental  difficulties  in  focusing  and  beam  deflection  and  would  have 


2 


required  modifications  to  the  apparatus  in  its  present  form.  Use  of 
35 

isotopically  pure  Cl  would  have  been  made  once  it  were  ascertained  that 
the  measurements  were  feasible. 

Prior  to  the  beginning  of  our  experimental  effort,  theoretical 
calculations  of  the  Na  + Cl”  and  K+  + Cl”  reaction  cross  sections  were 
made  by  Dr.  Olson.  For  collision  energies  between  0.1  eV  and  20  eV , which 
would  be  covered  in  our  experiments,  the  calculations  gave  neutralization 
cross  sections  in  the  10  - 10  ^ cm^  range.  This  range  is  one  to  two 

orders  of  magnitude  lower  than  that  which  was  characteristic  of  our  earlier 
measurements  on  atmospheric  species.  For  adequate  sensitivity  in  measuring 

cross  sections  of  this  small  magnitude,  we  hoped  to  have  beam  currents 

“7  + ~8  - 

exceeding  10  A for  Na  and  5 x 10  A for  Cl  . We  used  a thermionic 

source  for  Na+  which  gave  a beam  source  of  0.5  p.A.  This  is  adequate  for 

our  purpose.  Our  duoplasmatron  source  for  the  negative  ions  curned  out  to 

be  unsuited  for  stable  operation  in  the  presence  of  significant  amounts  of 

F or  other  F-containing  gases.  We  were  not  able  to  obtain  F“  beams  greater 
-9 

than  about  3 x 10  A;  thus  cross  section  measurements  were  not  possible  with 
our  existing  negative  ion  source.  It  would  be  worthwhile  to  investigate 
the  performance  of  alternate  ion  sources,  including  the  colutron  source, 
for  producing  F ; however,  such  an  investigation  of  negative  ion  sources 
would  require  new  resources  not  available  under  this  contract. 


THEORETICAL 

Our  theoretical  investigation,  which  was  quite  successful,  is  fully 
set  forth  in  the  attached  preprint  by  Dr.  R.  E.  Olson,  which  will  be 
published  in  Combustion  and  Flame. 
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Menlo  Park,  California  94025 

ABSTRACT 

Quantum  close-coupling  techniques  employing  pseudo-potential  derived 
coupling  matrix  elements  between  reactant  and  product  states  have  been  used 
to  calculate  the  cross  sections  for  the  reactions 

Na  + Cl  - Na+  + Cl* 

and 

K + Cl  ~ K+  + Cl* 

The  calculations  were  performed  from  threshold  collision  energies  to  5.0  eV , 

o 

allowing  the  reaction  rates  to  be  determined  for  temperatures  T s 5000  K. 

Analytical  forms  for  the  reaction  rates  are  given.  Representative  reaction 

o - 14 

rates  at  2000  K for  the  ionization  reactions  are  k(Na+Cl)  = 4.3  x 10 

3 -13  3 

cm  /sec  and  k(K  + Cl)  - 1.0  x 10  cm  /sec,  while  for  the  ion-ion 

+ -9 

recombination  reactions  the  rates  are  found  to  be  k(Na  + Cl  ) * 3.7  x 10 
cm? /sec  and  k(K+  + Cl  ) » 8.5  x 10  ^ cm^/sec.  It  is  believed  that  the 
calculated  rates  are  accurate  to  within  a factor  of  two  or  three. 


INTRODUCTION 


There  is  considerable  interest  in  determining  the  electrical  properties 
of  turbulent  afterburning  exhaust  plumes  of  solid  propellant  rocket  motors. 
Modeling  calculations  have  been  employed  to  calculate  the  radar  interference 
effects  that  arise  when  communication  with  rockets  involves  transmission 
through  these  plumes.  In  a recent  paper,  Jensen  and  Pergament^-  have  shown 
that  for  systems  where  potassium  is  a major  rocket  exhaust  impurity,  above - 
equilibrium  electron  concentrations  are  principally  formed  by  the  chemi- 
ionization  reaction 

K + Cl  - K+  + Cl*  (1) 

followed  by  the  electron  detachment  reaction 

Cl"  + H - HC1  + e"  . (2) 

However,  one  of  the  major  uncertainties  in  the  above  modeling  study  is  the 
ionization  reaction  rate  for  K + Cl,  which  is  known  only  to  an  accuracy  of 
one  or  two  orders  of  magnitude.  There  is  a similar  interest  in  the  Na  + Cl 
ionization  reaction. 

In  this  paper,  we  address  ourselves  to  the  problem  of  calculating  the 
cross  sections  for  the  ionization  reaction 

Aik  + Cl  - Alk+  + Cl"  , (3)  ' 

and  the  recombination  reaction 

Alk+  + Cl"  - Aik  + Cl  , (4) 


2 


where  Aik  represents  Na  or  K.  The  reason  for  including  the  recombination 

2 

reaction  (4)  in  the  calculations  is  that  experimental  data  exist  for 
o 

1800  K which  can  be  used  to  test  the  theoretical  approach.  The  recombination 
cross  sections  are  obtained  along  with  the  ionization  cross  sections  in  the 
same  calculation. 

For  the  theoretical  cross  section  determinations,  we  have  employed  the 

3 

quantum  close-coupling  method  and  coupling  matrix  elements  obtained  from 

4 

pseudo-potential  calculations  to  obtain  the  transition  probabilities  for 
reaction.  The  cross  sections  are  then  derived  from  the  calculated  transition 
probabilities  and  the  reaction  rates  are  obtained  by  averaging  the  cross 
sections  over  a Maxwellian  distribution  of  collision  velocities.  When 
possible  we  compare  the  theoretical  reaction  rates  to  existing  experimental 
data . 

THEORETICAL  APPROACH 

The  calculation  of  the  cross  sections  for  Reactions  (3)  and  (4)  are 
extremely  amenable  to  theoretical  attack  since  we  need  only  consider  three 
reactant  and  product  channels.  These  are  the  Aik  + Cl  ionic  channel  and 

two  neutral  states  arising  from  the  interaction  of  the  ground  state  alkaline 

2 2 

atom  and  chlorine  in  its  two  lowest  fine  structure  states:  Pi  and  P . 

i ? 

Higher  lying  excited  states  of  the  neutral  alkali-chlorine  atom  lie  above 
the  asymptotic  energy  of  the  positive  ion  (alk+) -nega tive  ion  (Cl  ) pair 
and  hence  to  first  order  do  not  interact  with  the  attractive  coulomb 
potential  of  Alk+  + Cl  . The  collision  mechanism  is  primarily  a curve 


crossing  transition  induced  in  the  regions  of  the  curve  crossings  between 
the  long-range  attractive  coulomb  potential  (-1/R)  of  the  ion  pair  Alk+  + Cl 
and  the  covalent  potentials  (-C  /R^)  of  the  neutral  species.  At  the  large 

(i  D 

internuclear  distance  where  the  curve  crossing  occurs,  the  covalent  potentials 
of  the  neutral  species  are  essentially  constant  and  may  be  accurately 
replaced  by  their  asymptotic  dissociation  energies.  The  curve  crossing 
distances  for  the  two  systems  are  given  in  Table  I. 

The  quantum  close-coupling  technique  for  obtaining  the  transition 
probabilities  as  a function  of  impact  parameter  b and  collision  energy  E 

3.5 

has  been  described  many  times  previously.  ’ The  input  parameters  for  the 

close-coupling  calculations  are;  (a)  the  internal  energy  change  AE  for 

the  reaction;  we  have  used  the  spectroscopic  levels  from  volume  I of 
6 

Charlotte  Moore;  (b)  the  initial  and  final  state  potentials;  the  potentials 
described  in  the  previous  paragraph  have  been  modified  to  include  a repulsive 
wall  so  that  there  would  be  no  unphysical  trapping  of  the  particles  by  a 
purely  attractive  coulomb  potential;  and  (c)  the  coupling  between  the  two 
channels.  In  the  close-coupling  calculations  it  is  necessary  to  include  the 
R dependence  of  the  coupling  matrix  element.  We  have  therefore  used  the 
functional  form  described  previously^  with  a modification  of  the  pre-exponential 
coefficient  so  that  we  could  reproduce  the  coupling  at  the  crossing  point 
to  obtain  the  pseudo-potential  values  of  • For  the  Na_C^  systems, 

the  form  is  given  by 


4 


H^(R)  = 8.98  arR  exp(-1.064  R) 


(5a) 


while  for  K-Cl  we  used 

Hn(R)  = 7.89  ?R  exp  ( - 1 .017  R)  . (5b) 

In  Eqs.  (5),  the  internuclear  separation  R is  given  in  units  of  ?\  while 
the  energy  is  in  units  of  electron-volts.  Since  the  pseudo-potential 
calculations  did  not  incorporate  the  fine  structure  states  of  Cl,  it  is 

necessary  to  include  in  the  matrix  elements  a factor  or  which  is  (3)  for 

2 . i 2 

the  Aik  + Cl(  P->)  state  and  (3)  for  the  Aik  + Cl(  P ) state. 

* h 

The  output  of  the  close-coupling  calculation  is  the  transition 
probability  for  reaction  P(b)  vs  the  impact  parameter  b.  The  total  cross 
section  for  a given  collision  energy  is  then  given  by  the  integral 

CO 

Q = 2ttS  j'  db-b-P(b)  , (6) 

0“ 

where  S is  the  statistical  weight  factor  for  the  reaction.  For  the 
ionization  reactions,  we  have  employed  j-j  notation  for  the  molecular  states 

formed  during  the  collision  and  find  that  the  appropriate  statistical 

2 

weight  S to  be  used  for  ionization  from  the  Aik  + Cl(  P3)  state  is  1/8, 

a 

2 

while  from  the  Aik  + C 1 ( P.)  state  it  is  1/4.  For  the  reverse  reactions, 
that  of  ion-ion  recombination,  the  statistical  weight  factor  is  unity. 

We  may  average  the  cross  sections  over  a Maxwellian  distribution  of 
collision  velocities  in  order  to  calculate  the  reaction  rates: 
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00 


(7) 
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(2Tru,kT) 2 


J dE-E-Q(E) • exp ( -E/kT) 
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Moreover,  since  the  recombination  cross  sections  at  Low  collision  energies 

g 

are  proportional  to  1/E,  a formula  for  the  Low  temperature  reaction  rate 
may  be  derived  from  Eq . (7)  that  varies  as  T 


k 


E'Q(E) J 
k 

(^nkT) 


(8) 


It  is  also  possible  to  calculate  the  cross  sections  for  Reactions  (3) 

and  (4)  by  using  the  Landau-Zener  method.  However,  this  method  cannot  be 

expected  to  yield  cross  sections  with  an  accuracy  greater  than  a factor  of 

9,10 

two  even  when  the  coupling  terms  are  known  exactly.  Furthermore,  for 

a curve  crossing  at  a large  internuclear  separation,  such  as  in  the  K-Cl 
system,  the  Landau-Zener  method  may  greatly  underestimate  the  cross  section 
due  to  the  neglect  of  the  coupling  over  an  extended  region  near  the  curve 
crossing  where  the  potential  energy  curves  are  still  close  together. 


RESULTS  AND  DISCUSSION 

The  calculated  ionization  cross  sections  for  reactions  (3)  are  given 
in  Table  II.  At  collision  energies  just  above  threshold,  the  cross  sections 
are  almost  constant.  This  behavior  is  expected  since  the  collision  velocity 
which  induces  the  transitions  and  determines  the  transition  probabilities 
is  varying  only  slightly  in  the  threshold  region,  and  because  the  maximum 
impact  parameter  collision  that  contributes  to  the  cross  section,  Eq . (6), 
is  also  approximately  constant.  We  should  reiterate  that  the  ionization 
cross  sections  given  in  Table  II  include  the  1/8  statistical  weight  factor 


6 


for  Aik  + Cl(  P3)  collisions  and  a 1/4  statistical  weight  factor  for 
2 

2 

Aik  + C 1 ( P, ) collisions. 

5 

The  ion-ion  recombination  cross  sections  for  Reactions  (4)  are 
presented  in  Table  III.  As  expected,^  the  low  energy  cross  sections,  which 
are  dominated  by  the  strongly  attractive  coulomb  potential,  decrease  as  1/E 
yielding  a low  temperature  reaction  rate,  Eq . (8),  that  is  proportional 
co  T-\ 

The  ionization  and  recombination  cross  sections  given  in  Tables  II 

and  III  may  be  used  to  calculate  the  reaction  rates  by  employing  Eq . (7). 

For  the  ionization  reactions,  two  rates  are  obtained  for  each  system  that 

2 2 

correspond  to  ionization  from  the  Aik  + Cl(  Pj_)  and  Aik  + Cl(  P ) fine 

structure  states.  For  low  temperature  reactions,  T ^ 1000°K,  the  reactant 

2 

Cl  will  be  preferentially  in  its  ground  P3  state  so  that  the  rate  given 
2 

for  Aik  + Cl(  P3)  will  apply.  However,  at  higher  temperatures,  the  fine 
T 

structure  states  of  Cl  will  most  likely  be  populated  statistically  with 
2 2 

5 in  the  Cl(  P3)  state  and  5 in  the  Cl(  P ) state, so  that  a weighted 
i h 

combination  of  the  rates  given  should  be  used. 

For  K + Cl  ionization,  we  find  the  following  reaction  rates: 

2 -14  v 

k[K  + Cl(  P|)]  - 2.2x10  T^ ( 1 + 8400/T)  exp  (-8400/T)  (9a) 


k[  K + Cl  ( 2P,  ) J =*  2.6  xl0'K  T^(l  + 7140 /T)  exp  (-7140/T) 
Likewise,  for  ionization  of  Na  + Cl,  the  rates  are  found  to  be 


and 


2 -13  if 

k[Na+Cl(  PO)  - 6.2x10  T (1  + 17700/T)  exp  (-17700/T)  (10a) 

i 


k[  Na  + Cl  ( ~ P ) ] = 4.3xlo'13  (1  + 16400/T)  exp  (-16400/T)  . (10b) 

3 , 

In  Eqs.  (9)  and  (10),  the  reaction  rates  k are  in  units  of  cm  /sec,  while  the 

o 

temperature  T is  in  K.  and  the  exponential  constants  are  the  thermodynamic 
endoergicities . The  temperature  range  of  validity  for  the  rates  is  T C 5000°5C. 
The  ion-ion  recombination  rates  are  calculated  to  be 


k(K+  + Cl~)  = 3.1  xl0~9  t"^  (1  + 2.6xl0‘6  T2) 


(lla) 


and 


■ • 7 • ir  7 — 

k(Na  + Cl  ) =»  1.6  x 10  T ^(1  + 4. 7 x 10  T2)  . 


(lib) 


The  units  are  the  same  as  those  given  for  the  ionization  rates.  In 
Eqs.  (11),  the  rates  are  predominately  a function  of  T with  a small 
correction  factor  for  the  higher  temperatures. 

We  may  compare  the  calculated  reaction  rates  to  two  sets  of  experimental 


measurements.  For  ion-ion  recombination,  Hayhurst  and  Sugden  have  indirectly 

o + 

measured  the  reaction  rates  for  Reactions  (4)  at  1800  K.  For  Na  + Cl 

-10  3 

recombination,  the  experimentally  derived  rate  is  6 x 10  cm  /sec,  while 

-9  3 + - 

the  theoretical  value  is  3.9x10  cm  /sec.  For  K + Cl  recombination 

-10  3 

there  is  a smaller  difference:  the  experimental  value  is  1.6x10  cm  /sec 

-11  3 

while  the  theoretical  rate  is  8.8x10  cm  /sec.  The  uncertainties  in  the 
experimental  rates  have  been  estimated  by  Jensen  and  Jones33  to  be 


approximately  a factor  of  5.  Moreover,  since  we  expect  the  theoretical 


values  to  be  accurate  to  only  a factor  of  two  or  three,  we  can  conclude 
that  the  theoretical  recombination  rates  are  in  reasonable  agreement  with 
experiment . 

No  direct  measurement  of  the  ionization  reaction  rates  exists;  however, 

a rough  value  for  K + Cl  ionization  has  been  inferred  from  the  modeling  of 

o 

the  transverse  radar  attenuation  in  a rocket  plume.  At  T s 1600  K,  Jensen 

1 -13  3 

and  Pergament  have  used  a value  of  k(K  + Cl)  = 1 . 3 x 10  cm  /sec  and 

conclude  that  this  value  is  too  large  by  approximately  an  order  of  magnitude. 

o 

In  comparison,  the  theoretical  value  at  1600  K for  a statistical  population 

-14  3 

of  the  Cl  fine  structure  states  is  4.1x10  cm  /sec.  Hence,  the  theoretics 

rate  is  smaller  than  the  original  K + Cl  ionization  rate  used  by  Jensen  and 

Pergament,  and  is  in  the  direction  to  more  accurately  reproduce  the  radar 

attenuation  data.  However,  our  rate  is  still  a factor  of  3.2  greater  than 

their  best  estimate  for  this  reaction  rate. 

We  have  also  used  a statistical  combination  of  the  ionization  reaction 

rates  along  with  the  ion-ion  recombination  rates  to  estimate  the  equilibrium 

constant,  k , where 
eq 


k ■ k / k 

eq  recomb  ionization 


(12) 


The  theoretical  equilibrium  constant  for  K-Cl  in  the  temperature  range 
1000°K  > T 4 3000°K  is  found  to  be 


theory 

k 1 (K  - Cl)  - 15.9  exp  (8180/T) 
eq 


(13) 
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This  expression  may  be  compared  to  that  given  in  the  JANAF  tables: 


IANAF 

k (K - Cl)  - 11.4  exp  (8200/T)  . (14) 

eq 

The  theoretical  equilibrium  constant  is  approximately  357.  larger  than  the 
form  given  in  the  JANAF  tables.  Similarly,  for  the  Na  -Cl  reactions  we 
obtain 

^theory  (Na_cl)  , 15>6  exp  (17200/T)  , (15) 

eq 

while  this  expression  may  be  compared  to 
JANAF 

k (Na  - Cl)  = 11.3  exp  (17480/T)  . (16) 

eq 

o 

In  the  temperature  range  of  1000  to  3000  K,  equilibrium  values  obtained 
from  Eqs . (15)  and  (16)  agree  to  within  25%;  In  both  Eqs.  (13)  and  (15),  the 
constants  were  obtained  from  a curve  fit  to  the  temperature  dependence  of  the 
theoretical  equilibrium  values. 

9 

We  have  also  used  the  Landau-Zener  (LZ)  method  to  estimate  the  ioniza- 
tion cross  sections  for  the  Na  + Cl  and  K + Cl  reactions.  For  Na  + Cl,  the 
LZ  cross  sections  are  in  reasonable  agreement  and  are  ~ 407.  lower  than  the 
close-coupling  results.  However,  for  K + Cl,  the  LZ  cross  sections  are  several 
orders  of  magnitude  below  the  close-coupling  results.  The  reason  is  that  the 
LZ  method  only  allows  transitions  at  the  crossing  point.  Thus,  for  a case 
with  a distant  crossing  such  as  K + Cl,  the  transition  probabilities  are 
greatly  underestimated  since  strong  coupling  can  occur  over  an  extended  region 
on  either  side  of  the  curve  crossing  where  the  potential  energy  curves  lie 

13 

close  together.  A similar  conclusion  has  been  arrived  at  by  Weiner  et  al. 
in  the  analysis  of  their  experimental  data  on  the  Na+  + 0 reaction. 


The  calculated  transition  probabilities,  P(b),  have  a regular  oscillatory 


structure  as  a function  of  impact  parameter  which  is  associated  with  a curve 
14 

crossing  process.  For  both  the  Na-Cl  and  K-Cl  systems,  the  magnitude  of 
P(b)  are  largest  at  small  impact  parameters  and  rapidly  decrease  to  zero  for 
impact  parameters  whose  classical  turning  points  are  outside  the  crossing 
distance.  We  can  obtain  an  estimate  of  the  average  transition  probability 
for  a specific  collision  energy  by  using  the  ionization  cross  sections  given 
in  Table  II  and  the  formula 

' qiON  + •’"x2  ' <"> 

If  we  are  concerned  with  the  transition  probability  for  scattering  only  on 

the  ground  state,  then  the  P(b)  calculated  from  Eq.  (17)  must  be  divided 

2 

by  the  statistical  weight  factor  S which  is  1/8  for  A1X+  Cl(  P.)  collisions 

1 

and  1/4  for  A1K+  CI(  P ) collisions.  Representative  values  of  P(b)  for 
scattering  on  the  ground  state  are  1-5%  for  the  Na-Cl  system  and  0.01-0.1% 
for  the  K-Cl  system. 

Probably  of  greater  interest  is  an  estimate  of  the  branching  ratio 
between  neutral  and  ionic  products  in  the  collisional  dissociation  of  NaCl 
and  KC1  molecules.  For  this  problem  we  only  want  the  transition  probability, 
p,  for  a single  passage  of  the  crossing  points.  A good  approximation  is  to 
use  the  Landau-Zener  formula  to  solve  for  p: 

P - 2p(l-P)  . (18) 

For  the  Na-Cl  and  K-Cl  systems,  the  transition  probabilities  are  small  so  we 
may  set  p - P/2. 

It  is  more  difficult  to  make  an  estimate  of  the  branching  ratio  for  a 

0 

Boltzmann  distribution  of  gas  whose  mean  temperature  is  less  than  the 


r— %•  n.  ■ 


endoergicity  of  the  reaction.  We  may,  however,  use  the  rates  given  by 
Eqs.  (9)  and  (10)  and  divide  by  the  mean  velocity  to  obtain  a cross  section 
for  use  in  Eq . (17).  Such  a calculation  was  performed  for  Na-Cl  and  we 
find  the  theoretical  branching  ratio  of  neutral  to  ionic  products  is 

within  a factor  of  three  to  the  branching  ratio  measured  by  Ewing  et  al.^ 

* 

We  also  substantiate  their  observation  of  a large  temperature  dependence 

of  the  branching  ratio  for  NaCl  collisional  dissociation  and  find  the 

neutral  to  ionic  ratio  decreasing  by  a factor  of  50  between  640Q°K  and 
o 

3000  K.  For  the  KC1  system  we  also  concur  with  the  findings  of  Ewing 
et  al.^  that  the  collisional  dissociation  primarily  produces  ionic  products 
and  find  theoretical  branching  ratios  that  are  approximately  two  orders  of 
magnitude  smaller  than  for  NaCl. 
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SUMMARY 


The  theoretical  reaction  rates  given  by  Eqs . (9),  (10)  and  (11) 
should  be  accurate  to  a factor  of  two  or  three  since  the  potential  energy 
curves  and  coupling  matrix  elements  are  known  to  a sufficient  accuracy. 
Moreover,  we  would  not  expect  rotational  coupling  (we  have  only  considered 
radial  coupling)  to  be  important  in  inducing  transitions  since  the  colli- 
sion velocities  are  quite  low. 

This  work  was  originally  started  to  meet  the  need  for  more  accurate 
values  for  the  ionization  reaction  rates  for  K - Cl  and  Na  - Cl  collisions. 
It  is  hoped  that  the  reaction  rates  and  their  temperature  dependencies 
given  by  Eqs.  (9)  and  (10)  will  be  an  aid  to  future  rocket  plume  modeling 
s tudies . 

Very  recently,  theoretical  calculations  have  been  performed  by 
16 

Arora  et  al.  for  the  K - Cl  and  Na  - Cl  ionization  cross  sections.  High 

energy  approximations  were  used  to  estimate  these  cross  sections  for 

collision  energies  E ^ 500  eV.  However,  the  low  energy  velocity  dependence 

and  magnitude  of  the  cross  sections  predicted  by  these  authors  is  in 

disagreement  with  the  results  given  here.  Moreover,  extrapolation 

of  their  results  to  thermal  energies  would  yield  ionization  cross  sections 
-20  2 

on  the  order  of  10  cm  and  predict  reaction  rates  substantially  smaller 
than  the  experimental  determinations. 
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Table  I 


CURVE  CROSSING  PARAMETERS 


System 

R -X 

X 

4 E-eV 

K + C1(2P2) 
2 

19.9 

0.724 

K + Cl(2p,  ) 

23.4 

0.615 

Na  + C1(2P3) 
2 

9.5 

1.523 

Na  + C1(2P^) 

10.2 

1.414 

13 
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Table  II 


IONIZATION  CROSS  SECTIONS  VERSUS  COLLISION  ENERGY 


System 


K + Cl(2Pi) 
z 


K + C1(2P,  ) 
3 


Na  + Cl(2pO 
2 


Na  + CI( 


E - eV 

2 

Q . -cm 
ion 

0.734 

6 . 1 (-18) 

0.749 

6 . 1(-18) 

0.774 

5 . 6 ( - 18) 

0.824 

5 . 9 ( - 18) 

0.974 

6 . 9 ( - 18) 

1.224 

8 . 3 ( -18) 

1.724 

1.1(-17) 

3.224 

1 -5 (-17) 

5.724 

1 • 8 ( - 1 7 ) 

0.625 

7 . 8 ( - 18) 

0.640 

9 . 1 ( - 18) 

0.665 

8 .0 ( -18) 

0.715 

7.5(-18) 

0.865 

7 . 9 ( - 18) 

1.115 

9 .5 (-18) 

1.615 

l.K-17) 

3.115 

1 . 4 ( —17) 

5.615 

1 • 6 ( — 1 7) 

1.533 

1.4(-16) 

1.548 

1 • 7 ( -16) 

1.573 

1 . 7 ( - 16) 

1.623 

1.4(-16) 

1.773 

1 .6 ( -16) 

2.023 

1 .6 (-16) 

2.523 

1 . 2 ( — 16 ) 

4.023 

1.1(-16) 

6.523 

9.9(-17) 

1.424 

9. 7 (-17) 

1.439 

7 - 3 ( — 17) 

1.464 

6 .6( -17) 

1.514 

7 . 3(— 17) 

1.664 

7 . 2 ( — 1 7) 

1.914 

6 . 1( -17) 

2.414 

6 . 5 ( - 17) 

3.914 

5 . 8 ( - 1 7) 

6.414 

5 . 7 ( — 1 7) 

*6 . 1 (-18) 


means  6 . 1 
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Table  III 


ION-ION  RECOMBINATION  CROSS  SECTIONS  VERSUS  COLLISION  ENERGY 


System  E-eV 


K + Cl  0.010 

0.025 
0.050 
0.10 
0.25 
0.50 
1.0 
2.5 

5.0 

Na+  + Cl"  0.010 

0.025 
0.050 
0.10 
0.25 
0.50 

1.0 
2.5 
5.0 

★ 

5.5(-15)  means  5.5  x 


, 2 
! -cm 

recomb 

2 

EQ-eV  cm4’ 

5 ■ 5 ( - 15) * 

5 . 5 ( - 1 7) 

2 . 2 ( - 15) 

5 - 6 ( - 1 7) 

l-K-15) 

5 • 6 ( - 1 7) 

6 .0 (-16) 

6 .0  (-17) 

3 . 2 ( - 16) 

8 . 1 (- 17) 

2 . 5 ( - 16 ) 

1 .21-16) 

2. 2(-16) 

2. 2 ( - 16 ) 

2. 2 ( - 1 6 ) 

5 .6 ( -16) 

2.4(-16) 

1 • 2 ( - 15) 

2 . 5 ( - 13) 

2.5(-15) 

9 . 9 ( - 14) 

2 . 5 ( - 15) 

4.9(-14) 

2.5 (-15) 

2 . 6 ( — 14) 

2 . 6 ( - 15) 

1.1(-14) 

2 • 7 ( - 15) 

6 .0  ( - 15) 

3 .0 (-15) 

3- 1 ( - 15 ) 

3. 1 ( - 15) 

1 . 8 ( - 15) 

4 . 6 ( - 15) 

1 • 3 ( — 15 ) 

6 .6  (- 15) 

-15 
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